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We investigate the spectrum and wave functions of qq' bound states for heavy fourth generation quarks (q') 
that have a very small mixing with the three observed generations of standard model quarks. Such bound states 
come with different color, spin and flavor quantum numbers. Since the fourth generation Yukawa coupling. A,/ , 
is large we include all perturbative corrections to the potential between the heavy quark and antiquark of order 
\gi Nc/IGtv'^ where A'^c is the number of colors, as well as relativistic corrections suppressed by iv/cf' . We 
find that the lightest fourth generation quark masses for which a bound state exists for color octet states. For 
the the color singlet states, which always have a bound state, we analyze the influence that the Higgs couplings 
have on the size and binding energy of the bound states. 



I. INTRODUCTION 

Among the mysteries of nature is the number of genera- 
tions. We observe three generations, however there could be 
a fourth generation if the masses of the quarks and leptons 
are beyond our present experimental reach. Data from the 
Tevatron (under some circumstances) restricts the masses of 
the t' and b' quarks in a fourth generation to begreater than 
about 350 GcV. For some recent studies see yj]. A strong 
constraint on the masses of fourth generation quarks comes 
from precision electroweak physics. Heavy fourth generation 
quarks contribute to the 5* parameter and to the p parameter. 
Their large contribution to the S parameter rules out a fourth 
generation with degenerate t' and h' quarks. However they 
also contribute to the p parameter and Kribs et. al. showed 
that an acceptable combined fit to precision electroweak data 
can be achieved, for example, with a mass splitting of about 
50 GeV between fourth generation quarks in the mass range 
350 - 700 GeV H. See also the earUer work in fl. Large 
splitting may also be possible j|4|]. For a more recent discus- 
sion of electroweak fits see 

In addition, there is a "unitarity upper limit" on fourth gen- 
eration quark mass of about 500 GeV IMS- This does not, 
however, necessarily forbid heavier quark masses. Rather it 
indicates that higher order perturbative coiTections become 
important at this mass ifioll . Dynamical considerations rather 
than unitarity give an upper bound of about 3 TeV. This is 



similar to the upper bound on the Higgs scalar mass 111 ll 112 1 



A heavy fourth generation can destabilize electroweak sym- 
metry breaking (see lIlSll for a review). According to the re- 
cent work in ll4ll . if there is no new physics (apart from the 
fourth generation) below a TcV, the Higgs mass should be 
roughly equal to or larger than fourth generation quark mass 
in order to avoid the instability. Of course there could be new 
particles beyond the fourth generation fermions below a TcV 
that get a large part of their masses from electroweak symme- 
try breaking. For example, scalars S that have a term in the 
scalar potential gS^SH^H get a contribution to the squares 
of their masses equal to [v ~ 246 GeV) and such in- 

teractions could help stabilize the Higgs potential. 

It is easy to imagine simple physical mechanisms that sup- 
press the mass mixing between the heavy fourth generation 



quarks and the three generations of standard model quarks. 
For example, the fourth generation quarks and leptons could 
have a different value for B~L than the standard three genera- 
tions. (Here B and L are baryon number and lepton number.) 
If B — L violation is small then the mixing between fourth 
generation quarks and the standard three generations is sup- 
pressed. For example, a fourth generation of quarks and lep- 
tons with both baryon and lepton number minus three times 
those of the ordinary three generations of quarks and leptons 
can cancel the baryon and lepton number anomalies allowing 
those symmetries to be gauged lITsll . 

Heavy fourth generation quarks feel a strong attractive 
force from Higgs exchange in both the q'q' and q'q' channels 
that gives rise to bound states lfl6ll . If the fourth generation 
quarks have a very small mixing with the ordinary quarks, 
they can be long enough lived that bound q'q' states decay 
through q'q' annihilation and not via q' decay to a lower gen- 
eration quark and a W boson. In this case the production of 
these bound states at the LHC may have important experi- 
mental consequences. Furthermore the n'q' bound states may 
be long very lived. References lll6l - ll9ll discuss some other 
interesting possible physical consequences of a heavy fourth 
generation. 

In this paper we focus on the physics of the q'q' states. Here 
we explore the role of perturbative coiTections suppressed by 
NcX'^, /IGtt^, and as (here Nc is the number of color and ag is 
strong coupling constant), as well as relativistic corrections on 
the the spectrum and wave functions of the q'q' bound states. 
We find that the perturbative and relativistic corrections have 
a significant impact on wave functions and spectrum of q'q' 
bound states. 

The q'q' bound states can be in a color singlet or color octet 
configuration. For the color octet states we find the lightest 
fourth generation quark masses for which a bound state ex- 
ists. In any color singlet configuration there is always a bound 
state. Therefore we discuss the impact of the Higgs couplings 
to the heavy quarks on the shape of the wave functions for the 
bound states and the bound state binding energies. In the nu- 
merical analysis, we sometimes show results for values of niqi 
that are below the experimental limit of 350 GeV or above 
500 GeV where we expect perturbation theory to be of lim- 
ited use. Our excursion in to these regimes is for pedagogical 
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FIG. 1 : Diagrams of Higgs self energy (left) and W boson vacuum 
polarization (right). 



reasons and does not mean we dismiss the constraints from 
experiment or the limitations imposed by perturbativity. 

II. HAMILTONIAN 

Since precision electroweak physics favors a small value for 
|(mf,/ — mti)/{mi,i + mt') \ (here m^/ and mf,/ are masses of 
t' and b', respectively), we work in the limit where the heavy 
fourth generation quark masses are equal, i.e., nif = irii,/ ~ 
iriqi . It is straightforward to add in the effects of the differ- 
ence between the heavy fourth generation quark masses. (We 
discuss the impact of a fourth generation quark mass splitting 
at the end of this paper) Then the leading order Hamiltonian 
for heavy quark bound states from Higgs scalar exchange is 



(\/2G 



(1) 



where p and x (r = |x|) are momentum and relative coordi- 
nate in the center of mass frame, and to^ is the Higgs scalar 
mass. In momentum space the leading potential from tree- 
level Higgs exchange is 



V{p) 



V2G 



(2) 



Here we have expressed the heavy quark Yukawa coupling, 
A„' , in terms of the Fermi constant G p and the heavy fourth 



J 



Here the functions L'^ and are defined by 



generation quark mass. Using the variational method based on 
a trial wave function ^I* oc e"*^/" llT6ll and taking, for example, 
mh = 130 GeV, this Hamiltonian has an 5'-wave bound state 
forforTOq' > 583 GeV. In this section, we compute perturba- 
tive corrections and relativistic corrections to the Hamiltonian. 
We also give the QCD potential at order of as for color singlet 
and color octet configurations. 



A. Perturbative Corrections to the Potential 
Enhanced by the Number of Colors 

Here we include perturbative corrections to the potential of 
q'q' state of order N c\^, / IQ'k'^ ■ These arise from the heavy 
fourth generation quark contribution to the Higgs boson self- 
energy Yj^{p'^). They will give a correction to the leading or- 
der potential. For pedagogical reasons we also include terms 
proportional to the top quark Yukawa squared but set the other 
quark masses to zero. (In the numerical results we will see the 
top quark contribution is negligible.) Expressing the Yukawa 
coupling squared in terms of the quark mass squared and Gp, 
all the perturbative corrections enhanced by a factor of Nc 
come from the Higgs scalar self-energy and the W boson vac- 
uum polarization. 

Let us consider the Higgs propagator It is determined by 
the one-loop calculation of quark loop diagram (see left on 
Fig.ffl), 



Dh ip ) - — =h : 



(3) 



with 



q=t',b',t 

I 2 2-,dLl:{p^) 



dp"^ 



(4) 
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2o3 



Ll{p^ 
L\{p' 



4m^ < 

2^253 tan^i ( i ) < p2 < 4m2 



(5) 



where 



1 - 



4m2 



b = 



/ 4m2 77^2 

1 and Xn = — T 



(6) 



The derivative of (p^) evaluated at p^ 



dp^ 



1-Axq + 2y/Axg - 1(1 + 2xg)tan-\l/y/4:Xq - 1). 



(7) 
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Finally, using MS subtraction (with the number of space-time dimensions n = 4 



e), S'' ; 



IS given as 



5'^ 



E 

q=t' ,b',t 



167r2 



log 



-2 



(8) 



in a power series in, — m\. 



Expanding the factor S [p 
is clear from its definition in Eq. (|4) that it first contributes at 
order (p^ 



2\2 



4) 

In the numerator of the Higgs propagator the factor 6'^ is 
divergent. This divergence is cancelled in the potential if we 
express the fourth generation quark Yukawa couplings A^, in 
terms of G fiti^, . The resulting correction to the Fourier trans- 
form of the potential is 



Vport(p) 



(9) 



FIG. 2: t-channel Higgs exchange. 



Here the perturbative corrections in the denominator of the 
Higgs propagator can be negligible. This arises because of a 
cancellation between the three terms in Eq. (|4]i. Here we have 
used the expansions. 



? 9' 



L<{ml) = ml 



dL'kip') 



dp^ 
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(10) 



(11) 



(12) 



Even though expansions we used in Eqs. (fT0l)-(fT2t are appli- 
cable for rrih <C 2mqi, we have checked numerically that 
the denominator of the Eq. (|9]l is a good approximation when 
nih ~ . On the other hand, the factor 6 in the numerator 
of the potential is given by 



(13) 



Here n^^(p2) is the transverse part of W boson vac- 



uum polarization Ilww,fj,iyip 



defined by Uww,f,,y{p^) = 



9iii^^wwiP^)~^' ■ ■ ■ nvi'iv(O) there are two contributions 



(see right on Fig. [Til, 11^ 



ni;.v^,.(o)+n^ 



these are 
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(0) and 



wwg'iO)/Mi 



167r2 

Using the above results, one obtains 
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FIG. 3: s- and f-channel neutral fictitious scalar exchange. 



Here we take Nc = 3. The divergences in 5^ and n^^(O) 
canceled. Eqs. ^ and (fT6l i are the main results of this section. 



B. Relativistic Corrections 

Relativistic corrections to the potential come from expand- 
ing the spinors in the Higgs scalar exchange diagram and from 
including the contributions from longitudinal gauge bosons 
and the "fictitious scalars" in gauge. We choose ^ = 1 so 
that the longitudinal gauge boson contribution vanishes. The 
coiTections appropriate for the ground 5-wave bound states 
are given here. For discussions of how the relativistic correc- 
tions to the potential are derived from Feynman diagrams see 
Refs. SIM]. 

Expanding the spinors for the ^-channel Higgs exchange di- 
agram (shown in Fig.|2]i gives the relativistic correction to the 
potential. 



K-ol. Higgs (p) = - 



V2Gi 



(17) 



Neutral fourth generation bound states can exist in the fla- 
vor states Pt' and b'b', can be in color singles (1) and octets 
(8), and furthermore they can have zero and one spins. Since 
we are working in the limit mj/ = nib' , it is convenient to de- 
compose the flavor structure into heavy quark isospin / = 
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p+ 



FIG. 4: f-channel charged fictitious scalar exchange. 



find that, 



MS), 



(20) 



w 



where Mw is boson mass. Here plus and minus sign cor- 
respond to / = 1 and channels, respectively. 

Finally there is the usual relativistic correction to the kinetic 
energy. 



p2 / 

-'ici ^ — T 

■niqi \ Ain't, 



(21) 



FIG. 5: Gluon exchange. 

(;.e., {t'f + b'b')/V2) and / = 1 (i.e., {t'f - b'b')/V2). So 
far the contributions to the Hamiltonian have not depended 
on the bound states heavy quark isospin, color and spin quan- 
tum numbers. However, the contributions we consider now 
do depend on these quantum numbers. We therefore attach 
the superscript, (C = color, I = heavy quark isospin, S = 
heavy quark spin) to the potential. Since these eight states 
characterized by (C, /, S) do not necessarily form bound 
states, hereafter we call them "channels". 

Exchange of the neutral " fictitious scalar" (which we call 
P°) in the ^-channel (left on Fig. O gives spin-dependent po- 
tential, but is independent of the color and flavor We find 
that. 



Pf ,t-channel 



(P) 



Ml 



^(S), (18) 



where Mz is Z boson mass, and for spin one, 51(1) — —1/3, 
and for spin zero, 17(0) = 1. This contribution is attractive 
in the spin one channel and repulsive in the spin zero channel. 
The s-channel P° exchange (left on Fig. [3]), on the other hand, 
gives a repulsive potential which only occurs in the (1, 1, 0) 
channel. 



f>(ia,o) , ^ 

^P".,s-channelvP/ 



1 



Mll{Aml,y 



It is enhanced by a factor of = 3, compared with other 
relativistic corrections to the potential. 

The final relativistic correction to the potential comes from 
t-channel exchange of the charged fictitious scalar P+, which 
is depicted in Fig.|4] It is independent of color but depends on 
spin and on flavor since it mixes the t't' and b'b' channel. We 



C. QCD Potential 

There are also contributions to the potential from one-gluon 
exchange (Fig. |5]). They are attractive in the color singlet 
channel and repulsive in the color octet channel but are spin 
and flavor independent, and given as 



(1) 



QCD 



Vk%{r) 



(r) = 



"QCD 



6 



(22) 
(23) 



There are always bound states in the color singlet channel 
because the strong interactions confine. In the octet channel 
there are no bound states without the Yukawa potential from 
Higgs exchange. In our numerical work we evaluate at the 
Z boson mass, as{Mz) = 0.118. 

III. NUMERICAL RESULTS 

In this section we discuss the ground state S'-wave states in 
the various color, heavy flavor isospin, and spin channels. The 
Hamiltonian for this system is 



(24) 



where 



H^^^ — Ti-cl + ^pcrt + ^-ol. Higgs + t.(,,j3j,„gi 

+ ^P0, s-channel + ^P+ + VqcD • (25) 

We use the variational method, minimizing E[a] = 
{ij)\H\tp) / for trial wave functions oc e^''/". In order 
for the V jc expansion to make sense, we restrict our analysis 
(19) to wave functions that give an expectation value for p'^/m^, 
that is smaller than 1/3. This ensures that higher order terms 
in the vjc expansion, which we have neglected, are not im- 
portant. This means that. 



q'- 



(26) 



Before discussing the numerical results, we give the for- 
mula for E[a] in each channel. The expectation values of the 
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kinetic energy and the potential from the Higgs exchange and contribution for color singlet/octet and isospin zero/one chan- 
t-channel neutral fictitious scalar exchange give a common nels. These are given by. 



1 



IJlq/a 



+ S — to^/4to^, 
(2 + amh)^ 



n{s) 



7 9 



(2 + aMzf 



(27) 



The first term comes from the kinetic energy, while first and second terms in the second parentheses are from t-channel Higgs 
exchange, including the perturbative correction to the Higgs propagator The rest is from neutral fictitious scalar exchange in 
f-channel. The s-channel neutral fictitious scalar exchange, on the other hand, gives a contribution only for the (1, 1, 0) channel, 
which is 



s-channel [^-1 



V2G 



1 



■na (1 - A/|/W,) mia 



(28) 



As we mentioned, this term always contributes as a positive (repulsive) term in total energy, and it is enhanced by color factor 
Nc = 3. Charged fictitious scalar exchange gives 



E^P%] 



V2G 



pm^ 



na 



1 



2TO^,a^ 



(2 + aMw)' 



(29) 



Finally, the contribution from one gluon exchange is 



•^QCD ["] 



3a 



6a 



(30) 



With all terms we have given above, the variational energy in 
each channel is obtained, 

i?(^."-^)[a] = i?£), + 4^?+<)^, (31) 

i?(^'i'0)H = E('J^-E^^l+E^^Vf+E^^l^, (32) 

i?(^'i'i)[a] = Eill-E^^l+El^^^, (33) 
for color singlet state, (1, /, S), and 

£;(«-"^^)[a] = (34) 

E^'^'^'Ha] - (35) 

for color octet state, (8, /, S). These results are summarized 
in the Appendix. 

We compute the variational energy [a] in each 

channel and study the properties of the bound states. In the 
color singlet channels, there always exists a bound state. For 
small enough niqr the state is very close to the familiar QCD 
"onium" states. However as iriq' increases the parts of the po- 
tential proportional to m^, become more important. We find 
the value of a (in the parameter region given by Eq. (EFt ) 
which gives minimum binding energy for fixed iriqi. (We 
denote call it gq.) It is compared with Bohr radius of pure 
QCD potential, aqcD = 2as/3mq>. We begin by taking 
mil = 130 GcV. Later we redo the analysis for the case 
rrih — ruqi which may provide more realistic values of the 



Higgs mass given the constraints from stability of the Higgs 
potential. The results (for rrih = 130GcV) are shown in the 
upper panel of Fig. |6] From this figure, it can be seen that the 
size of the singlet bound states are not close to a QCD-like 
bound state when m.qr > 400 GcV, except for the (1, 1, 0) 
channel. The sharp break in behavior as rjiqi increases is due 
to the limit we impose on the value of gq {i.e., it is greater 
than or equal to \/3/mq>), which ensures that relativistic cor- 
rections are not too large . In the (1, 1, 0) channel, the con- 
tribution from the repulsive s-channel exchange potential 
is so large that the bound state has uq > aqcD for a range of 
masses. In the lower panel of Fig. |6] we plot the variational 
binding energy computed at a = ao for each color singlet 
channel. We find binding energies of O((10 — 100) GcV) for 
mq> - 400 - 500 GcV. 

For the color octet channels, on the other hand, bound states 
do not exist if the heavy fourth generation quark is too light 
and of course the Higgs Yukawa couplings always play a cru- 
cial role because the QCD potential is repulsive. In our nu- 
merical analysis, we find the lowest value of niq' for which 
the minimum of the variational energy E[a] (in the region, 
a > \/3/mqi) has a negative value. The results are summa- 
rized in Table U Note that the values of the fourth generation 
quark masses relevant here are not the ones in parenthesis. 
We find that the lower limit reduces to 440 — 570 GeV, com- 
pared to the one given by the leading order Hamiltonian (i.e., 
583 GcV). As in the color singlet channel, we plot the lowest 
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FIG. 6: (i) ao /oqcd as the function of quark mass in color singlet 
channels. Here we take rrih ~ 130 GeV and ao is the value for 
which E[a\ is minimized (and negative) for fixed quark mass, (ii) 
Variational binding energy of color singlet channels as the function 
of the heavy quark mass. Here we set a = ao- 



FIG. 7: Variational binding energy of color octet channels plotted 
as the function of the heavy quark mass. In the plot, we use irih = 
130 GeV and take a as the value which gives the lowest binding 
energy for fixed rrigi. Note that (8, 0, 1) and (8, 1, 0) channels give 
almost the same results. 
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(C, /, S) Lower limit of TUgi 

(8, 0, 0) 574 GeV (574 GeV) 

(8, 0, 1) 440 GeV (359 GeV) 

(8.1.0) 440 GeV (359 GeV) 

(8. 1. 1) 510 GeV (439 GeV) 



TABLE I: Lower limit of quark mass for which a bound state forms 
in the various color octet channels. We take = 130 GeV. The 
values in parentheses are given by using a > ypijraqi instead of 
Eq. lllll 
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variational binding energy for fixed m^/ in Fig. |7] This figure 
indicates that color octet bound states with binding energy of 
(9((10 - 100) GcV) exist when m,- ~ 450 - 550 GcV. The 
color octet states we found form color singlet hadrons by neu- 
tralizing their color charge at long distances with gluons and 
light quark-anti quark pairs. 

It is important to remember that when ao is at the end of 
the range given by Eq. (|26] |. the actual bound state may be 
relativistic and more deeply bound than the results presented 
in this section indicate. Such a situation occurs for the color 



FIG. 8: The same as Fig.|6l except for taking rriK — rrigi 

octet results, except for the (8, 0, 0) state, and in some of the 
color singlet channels at larger heavy quark masses. In order 
to see how our results are affected by the choice of region 
for a, we consider, for example, the case where expectation 
value of p^/to^, is less than 1/2, which corresponds to a > 
\/2/ ruqi . In Table|I] the lower limit on niqi for an octet bound 
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(C, I, S) 


Lower limit of m^/ 


(8,0,0) 


No bound state 


(8,0,1) 


534 GeV 


(8,1,0) 


534 GeV 


(8,1,1) 


696 GeV 



TABLE IL The same as Table. Ul except for taking rrih = 
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Introducing the notation, m± = rrit' ± TOts the energy eigen- 
values E± and mixing parameters B± are. 



E± 
B+ 



£;(C,0,S) £;(C,1,S) 




m_ 



m?., (38) 
(39) 
(40) 



FIG. 9: The same as Fig.l?] except for taking mu = 



respectively. When ao ^ aqcD, ti^ is larger in magnitude 
than (£'(C'0''S')-i<;(c,i,s)-)^2. (Here we are assuming \m- \ ^ 
50 GeV.) Then, the mixing parameter is not neglegible. On 
the other hand when ao is much smaller than aqcD, is not 
important and the mixing parameter is negligible. Then the 
states |±) are almost isospin eigenstates, i.e., |+) ~ |/ = 0) 
and |— )~|/=l).Ina more accurate evaluation, one should 
also take into account the correction m_ makes to £'('^'"'■5) 
and i?'^'-^'^'"^'. These corrections are suppressed by, (m_/7n+) 
and are expected to change the binding energies by a few to 
10%. 



state to exist is also derived using this region for a, instead of 
Eq. ( |26] | (see the values in parentheses). As it is shown, the 
limit becomes smaller by 10 — 20% (except for the (8, 0, 0) 
state). 

Finally, we give the numerical results in the case of mh = 
iriqi . In the same manner as we did when mh was fixed at 
130 GcV, ao/aqcD and the binding energy for color singlet 
channels are given in Fig.|6]and the lower bound on the fourth 
generation quark mass, and the binding energy for color octet 
channels are given in Table and Fig. |9] respectively. For 
singlet states, the bound state is certainly not QCD-like when 
■niq' > 400 (535) GeV in the (1,0, 1) ((1, 1, 1)) channels. 
On the other hand, octet channels do not form bound state 
unless niq' > 535 GeV, which is near the unitarity bound or 
equivalently strong coupling regime. 

We have assumed that mf = nib' throughout this paper, 
however it is straightforward to take into account the mass 
difference between the the heavy fourth generation quarks. In 
that case, the / = and 1 sates are no longer the energy 
eigenstates. Rather we denote the eigenstates by |+) and |— ). 
They are the following linear combinations of the heavy quark 
isospin eigenstates, 

1+) cx |/ = 0)+S+|/ = l), (36) 
I-) cx S_|/ = 0) + |/ = l). (37) 



IV. CONCLUDING REMARKS 

Heavy fourth generation quarks may have a long enough 
lifetime that it is sensible to consider their bound states. At the 
LHC heavy quark q'q' bound states will be produced by gluon 
fusion. Hence it is important to understand the properties of 
these states. In this paper we have determined the binding en- 
ergies and sizes of these states. For iriq/ > 400GeV, the Higgs 
Yukawa coupling plays a crucial role in the properties of these 
states and also relativistic and perturbative corrections are im- 
portant. In a future publication we hope to elucidate more 
of their properties, including production rates at the LHC and 
their decay branching ratios. 
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Appendix 



Here we give explicit formulas for the variational energy in each channel. 
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(1 + 5)a2 1 TO^ 

" 2(2 + am,,)2 + 4^ 8^ (2 + amnf 



Ml 



8to2, (2 + aMzf 4m2, (2 + aMwf 



3a 



a 



2\/2G 



(l + <5)a^ 
'2(2 + aTOh)2 



8m^ (2 + a?7ift)2 
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